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1 Introduction 



This is the first of a series of papers devoted to the study of hyperbolic geometric flow and its 
applications to geometry and physics. Hyperbolic geometric flow was first studied by Kong and 
Liu in [11) . To introduce such flow we were partially motivated by the Einstein equations in general 
relativity and the recent progress in the Hamilton's Ricci flow, and by the possibility of applying 
the powerful theory of hyperbolic partial differential equations to geometry. Hyperbolic geometric 
flow is a system of nonlinear evolution partial differential equations of second order, it is very 
natural to understand certain wave phenomena in nature as well as the geometry of manifolds, in 
particular, it describes the wave character of the metrics and curvatures of manifolds. We will see 
that the hyperbolic geometric flow carries many interesting features of both the Ricci flow as well 
as the Einstein equations. It has many promising applications to both subjects. 

The elliptic and parabolic partial differential equations have been successfully applied to dif- 
ferential geometry and physics. Typical examples are the Hamilton's Ricci flow and Schoen-Yau's 
solution of the positive mass conjecture. A natural and important question is if we can apply the 
well-developed theory of hyperbolic differential equations to solve problems in differential geome- 
try and theoretical physics. This series of papers is an attempt to apply the hyperbolic equation 
techniques to study some geometrical problems and physical problems. One has found interesting 
results in these directions, see for example |16j for the applications of the hyperbolic geometric 
flow equations to physics. Our results already show that the hyperbolic geometric flow is a natural 
and powerful tool to study some important problems arising from differential geometry and general 
relativity such as singularities, existence and regularity. In this paper we study the basic properties 
of the hyperbolic geometric flow such as the short-time existence, nonlinear stability and the wave 
feature of the curvatures. In the sequel we will study several fundamental problems, for example, 
long-time existence, formation of singularities, as well as the physical and geometrical applications. 

Let ^# be an n-dimensional complete Riemannian manifold with Riemannian metric gy, the 
Levi-Civita connection is given by the Christoffcl symbols 




{ 



dx % dx? dx l 




where g lJ is the inverse of gij . The Riemannian curvature tensors read 




dx l 




+ 




jp^iV Rijkl — dkpRiji- 



The Ricci tensor is the contraction 



Rik = ff* Rijkl 
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and the scalar curvature is 



R = g v R 



i j ■ 



The hyperbolic geometric flow under the consideration is the following evolution equation 



dt 2 



2R i3 (1.1) 



for a family of Riemannian metrics g%j{t) on ^£ . More general hyperbolic geometric flows were also 
introduced in [llj . A natural and fundamental problem is the short-time existence and uniqueness 
theorem of hyperbolic geometric flow (1.1). In the present paper, we prove the following short-time 
existence and uniqueness theorem, the nonlinear stability theorem for Euclidean space, and derive 
the corresponding wave equations for the curvatures. These results were announced in Kong and 
Liu [TT]- 

Theorem 1.1 Let , g^(x)) be a compact Riemannian manifold. Then there exists a constant 
h > such that the initial value problem 

-Qj?zgij(x : t) — 2i?y (a;, i) , 
,•(1,0) ^(1,0) = *&•(*), 

has a unique smooth solution gij(x,t) on j$ x [0, h], where kfj(x) is a symmetric tensor on .M . 

The main difficulty to prove this theorem is that, the hyperbolic geometric flow (1.1) is a system 
of nonlinear weakly-hyperbolic partial differential equations of second order. The degeneracy of 
the system is caused by the diffcomorphism group of ^# which acts as the gauge group of the 
hyperbolic geometric flow. Because the hyperbolic geometric flow (1.1) is only weakly hyperbolic, 
the short-time existence and uniqueness result on a compact manifold does not come from the 
standard PDEs theory directly. In order to prove the above short-time existence and uniqueness 
theorem, using the gauge fixing idea as in the Ricci flow, we can derive a system of nonlinear 
strictly- hyperbolic partial differential equations of second order, thus Theorem 1.1 comes from the 
standard PDEs theory. On the other hand, we can reduce the hyperbolic geometric flow (1.1) to 
a quasilinear symmetric hyperbolic system of first order, then using the Friedrich's theory [5] of 
symmetric hyperbolic system (more exactly, the quasilinear version [3]) we can also prove Theorem 
1.1. 

Noting an important result on nonlinear wave equations (see [10j). we will find its other inter- 
esting application to geometry, by applying this result to the wave equations of curvatures. In the 
present paper we first use it to prove the nonlinear stability of the flat solution of the hyperbolic 
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geometric flow defined on the Euclidean space with dimension larger than 4. More precisely, we 
have 

Theorem 1.2 The flat metric gij — dij of the Euclidean space M" with n > 5 is nonlinearly stable. 

See Section 4 for the precise definition of nonlinear stability which is very important in general 
relativity. The key point of the proof of this theorem is the global existence of classical solutions 
of the Cauchy problem for the nonlinear wave equations. 

Similar to Hamilton [5J, we derive the corresponding wave equations for the curvatures, for 
example, we have 

Theorem 1.3 Under the hyperbolic geometric flow (1-1), the curvature tensor satisfies the evolu- 
tion equation 

d 2 

-QpRijki — ARijki + 2 (Bijki — Bijik — Bujk + Bikji) 

—g Pq (RpjklRqi + RipklRqj + RijplRqk + RijkpRql) (1.2) 

+ 2 „ ( dT l dT % _ dT ldT%\ 

+ 9pq \ ot at dt dt )' 

where B^i = g pr g qs R p i q j R r ksi and A is the Laplacian with respect to the evolving metric. 

The wave equations for the Ricci and scalar curvatures are stated and proved in Section 5. This 
is similar to the Ricci flow equation, except that there are quadratic lower order terms involving 
the connection coefficients. It turns out that there is a very rich theory in nonlinear wave equations 
to deal with such terms, see |10) . 

From the above results we can already see that the hyperbolic geometric flow has many features 
of the Ricci flow, therefore many well-developed techniques in the Ricci flow may be applied to 
the study of this new kind of flow equations. On the other hand, the hyperbolic geometric flow 
can also be viewed as the leading terms of the vacuum Einstein equations. Since the hyperbolic 
geometric flow contains the major terms in the Einstein equations, it not only becomes simpler 
and more symmetric, but also possesses rich and beautiful geometric properties. See Section 6 for 
more detailed discussions on the relations between the hyperbolic geometric flow and the Einstein 
equations, and more generally its relations with other important problems in general relativity. 

The paper is organized as follows. In Section 2, using the gauge fixing idea as in the Ricci flow, 
we derive a system of nonlinear strictly-hyperbolic partial differential equations of second order. 
In Section 3 we reduce the hyperbolic geometric flow (1.1) to a quasilinear symmetric hyperbolic 
system of first order, and give the proof of Theorem 1.1. Section 4 is devoted to the nonlinear 
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stability of the hyperbolic geometric flow defined on the Euclidean space with dimension larger 
than 4. In Section 5, we derive the wave equations satisfied by the curvatures, and illustrate the 
wave character of the curvatures. Some discussions are given in Section 6. 



2 Strict hyperbolicity of hyperbolic geometric flow 

In this section we consider a modified system of evolution equations of the hyperbolic geometric 
flow, which is strictly hyperbolic so that we can get a solution for a short time by solving the 
corresponding Cauchy problem. The solution of the system (1.1) then comes from the solution of 
the modified equations. 

Let ^# be a compact n-dimensional manifold. We consider the hyperbolic geometric flow (1.1) 
on , that is, 

0^9ij{x,t) = -2R ij (x,t). (2.1) 

Suppose (jij(x, t) is a solution of the hyperbolic geometric flow (2.1), and ip t '■ *dt — ► ^ is a family 
of diffcomorphisms of Let 

gij(x,t) = i>tgij(x,t) 

be the pull-back metrics. We now want to find the evolution equations for the metrics g%j{x,t). 
Denote by y(x, t) = <pt(x) — {y 1 {x, t), y 2 {x, t), ■ ■ ■ ,y n (x,t)) in local coordinates. Then 

v^) = %%9^{y,t) (2-2) 



and 



d . . d 

dt 9ij{x ' t] = m 



, . dy a dy 13 



dy a dy 13 d , . . . , . d ( dy a dy 

^^dM y{x ^ t)+9 ^ t] dl (dx^dxl 



Furthermore, we have 



d 2 . dy a dy? d 2 g afj ( d 2 y<*\ dy^ 



+ £ A ( d ~X-) *«* + 2^ S d ~¥ (2-3) 



dx i dxi \ dt 2 J yap dx* V dt J dxi dt 



(dy a \ 






V dt j 


dxi 





On the other hand, 



+2 dxi g x j \ dt ) dt Jr2 dx l I ^ i ' ' 9af1 ' 



-(y{x,t),t) = -rr—^r + 



dt vyv ' n ' dy-y dt dt 



d 2 g a p ( , s , d 2 g a0 dy' 1 dy x d 2 g a pdy 1 d 2 g a(i dg afj d 2 y~> 

dt 2 w x > u > dyidyx dt dt dy^dt dt dt 2 dy dt 2 
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and 

It follows from (2.3) that 



dt 2 



-(x,t) 



2R af3 {y,t) 



dy a dyf 3 d 2 g a p dy a dyf 3 dy 1 dy* 
dx % dxi 



2 d 2 g a p dy a dy 13 dy 1 
dyfdt dx 1 dxi dt dx 



dyfdy x dx 1 dxi Qt 
d ( „ dy 13 d 2 y a 



9al3 



8x3 dt 2 



dt 



+ 



8x3 



9a(j 



Q y Q2 ya 



dx 1 dt 2 



dg a p dy a dy 13 



d 
dx 1 



dy' 1 dx 1 dx 3 Q x l \dx j;}>J1 J dxi \dx l 



d 



+2 



dx 1 



90i 

dg a /3 dy 1 



(dyf 3 , 

a 



dT\ dyP_ f dg af3 
dt J 8x3 \ dt ^ dyf dt 



90i 
dy 2 d 



9V 



dx 1 dx3 \ dt 



( dg a[j dy 1 dg a p 
\ dyf dt dt 



d 

' dx 1 



(dy a ^ 






\~dT; 


dxi 





Let us choose the normal coordinates {x l } around a fixed point p£ J' such that 
p. We next prove that, at p £ 



dgij 
dx k 



dg a p dy a dy 13 _ _d_ 
dx 1 



dyf 3 . 



90i 



dy 1 dx i dxi dx 1 \dxi 
The left hand side of (2.5) is 

dg a dy a dy 13 d V - 

dyf dx 1 dxi dx 1 { dx3 9Pl 



d (dy? „ 
dxl {dx^ 9 ^ 



_ _d_ (dy? . 
" dx3 {dx i9af3 



0, Vu7 = l, 



. n. 



(2.4) 
= at 

(2.5) 



d 

dr 



9mj 



dx m dx n \ dy a dy 13 



dy a J dyf 3 dx 1 dxi 



— gmj 



— gmj 



dy a dyf 3 J dx 1 dxi 
d (dx n 

d 

dx 1 V dyf 
d ( dx m \ dy a 
dx 1 



d 

dx 1 
dc 



Qmj 



dx" 
dy~i 



d 

dxi 



dx r 



dy 



{ dx n \ 


dy a 


dyf 3 ' 


{dtfj 


dx i 


dx3 _ 



dx r ' 



d f dx 



dyf 
d 

dy* V dy 



dy a 
dx m 



l dx3 \dyf 
d f dx 



- 9n 



dyf ydyf 3 
dy a d fdx m 

~dx 1+9mi dy7 3 {~dlf 



dy 13 
dxi 
dyf 



9m j 



d (dx" 



dx 1 \ dyf 
d (dx rn 



dx3 9mj dxi 



dyf 



- 9n 



d (dx m 
i dxi {~dyf 

d (dx m 
i dx~3 {~dyf 



0. 



So (2.5) holds. 
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By (2.4) and (2.5), we have 



dt 2 



-(x,t) 



. , d ( dx m d 2 y a 
. 2Rij{x ,t) + —{g mj — — 



d 

dx~J 



dy a d 2 y c 



dx m dt 2 



d 2 g afj dy a dy p dy 1 dy x t ^ dg a[j dy 1 dy a dy 13 



dy~<dy x dx l dxi dt dt dy^dt dt dx l dxi 



dx l V dt ) dxi \ dt + dy-y dt 



dx % dxi \ dt J \ dt 



lap 



+ 



dx l 



fdy a ^ 


A 


m 


Vdtj 


dxi 





Jg a p dy' 1 
dyt dt 

(2.6) 



We define y(x,t) — (ft(x) by the following initial value problem 

f «2 a a a o 



dt 2 dx k - 

y a {x,0) = x c 



d_ 

dt 



(2.7) 



y a (x,0) =yf(x) 



and define the vector filed 



1 jl L ji I , 



where T k - L and Y k t are the connection coefficients corresponding to the metrics gij (x, t) and g^ (x, 0) , 
respectively, yf(x) (a = 1, 2, • • • , n) are arbitrary C°° smooth functions on the manifold We 
get the following evolution equation for the pull-back metric 

<^"M - -2 R „(x, ( ) + V ( V J + V j V i+ 



+2-. 



dyidy x dx l dxi dt dt 
d 2 g al3 dy a dy 13 dy 1 d f dy^ \ dy^_ f dg afj dg a p dy 



+2 



dy^dt dx l dxi dt dx l \ dt J dxi \ dt dy~/ dt 
d_ ( dyP\ dy^_ ( dg afi dg al3 dy" 



dxi \ dt J dx l 



dt 



dyf dt 



(2.8) 



fdy a ^ 








dxi 





9ij(x,0) 



4 ^R^x, t) + ViVj + VjVi + F(Dy, D t D x y), 
d 

-gij(x,0) = fcO-(x), 



dt 



where 



Let 



Dy 



dy a dy a 
~~dT' ~dx 1 



X = 



D t D x y = 
dy a dy a d 2 y a 



d 2 y a 
dx l dt 



(a,i = l,2, ••• ,n). 



(a,i = 1,2, • • • ,n). 



dt dx % dx l dt / 

The nonlinear term F = F(X) = F(Dy, D t D x y) in (2.8) is smooth and F(X) = 0(\X\ 2 ) holds. 
Since 

pfc _ dy a dy 13 dx k ~ 7 ^ dx k d 2 y a 



dxi dx i dyf Q/3 dy a dxidx v 
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the initial value problem (2.7) can be written as 

d 2 y a = JX ( d 2 y a ° k d]f ™ V d]P_ 
dt 2 9 \dx3dx l ll dx k ^dx^dx i 

Q 

y a (x,0)=x a , —y a (x,0)=yf(x). 
At the same time, in the normal coordinates {x 1 }, 

■^M+w i+ ?,« - ^f&t-LW** ■ ds " 



(2.9) 



dx l \ dxi J dx k I \ 9a; 1 i9a; J eta' 



1,31,9 dx i \2 y \ dxi ^ dxP dx l 



+gikg j 



9 2 3fez # 9 5ii d 2 gij 



' dx l dxi dx l dx k dx^dx k dx k dx l 



2' [9a; z 9a;' 3 dx l dx q dx l dxi 
_1 Pg / 52 5p» , 5 2 si 9i 9%, 



2' [9a; 1 9X 1 ? dxidx? dx l dx^> 
g kl t ' J , + (lower order terms). 



(lower order terms) 



9a; fc 9:E i 

Thereby, the initial value problem (2.8) can be written as 



^fRz, = ^^1? + ^P*. D t D x y) + G(g, D x g), 
9 

gij(x,0) = g%{x), —gij(x,0) = ty(x), 
dgtj 

dx k 



(2.10) 



where 5 = D x p = ^-^ ) (i, j, = 1, 2, • • • ,n). Let 



dx k 

M= (s^ttt ) (i,3,k= 1,2,- •• ,ra). 



The nonlinear term G = G(£t) = G(g, -D^g) in (2.10) is smooth and quadratic with respect to D x g. 

We observe that both (2.9) and (2.10) are clearly strictly hyperbolic systems. Since the equa- 
tions (2.9) and (2.10) are strictly hyperbolic and the manifold ^0 is compact, it follows from the 
standard theory of hyperbolic equations (see [5], [5], [TO]) that the system united by (2.9) and 
(2.10) has a unique smooth solution for a short time. Thus, we have proved Theorem 1.1. 

3 Symmetrization of hyperbolic geometric flow — second 
proof of Theorem 1.1 

In this section we reduce the hyperbolic geometric flow (2.1) to a symmetric hyperbolic system. 
Then we use the theory of symmetric hyperbolic system to give another proof of Theorem 1.1. 
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Let be a compact 71-dimensional manifold and gij{x,t) is a hyperbolic geometric flow on 
^ . We denote the corresponding connection coefficients, the Riemannian curvature tensor and 
the Ricci curvature tensor by T k j, Rijki and Rik, respectively. 

We consider the space-time R x ^ with the Lorentzian metric 



ds 2 = —dt 2 + gij(x,t)dx i dx-' 



(3.1) 



and denote the corresponding connection coefficients, the Riemannian curvature tensor and the 
Ricci curvature tensor by T^, R a p-y\ and R a p, respectively. Here and hereafter, the Greek indices 
run from to n, Latin indices from 1 to n. The summation convention is employed. We also 
denote x° = t. 

By direct computations, 
1 dg io 



13 



pfc pO 

ij ' ij 



2 dt 



1 01 — 1 ZO — r,y Q. 1 1 Oi — U > 1 00 ~ U l I 



<9t 



00 



0. 



R 



R 



oka 



Then, 



dT% or* 



p£; pa pA? pa 



ijl i)r dxi ia jl ja a 
1 dg li dg ik , 1 u d 2 g lk 1 



R i3l + ~ A 9 



at at 



2 dt dt ' 2 9 



J n " n n 

dfi + 4 9 9 



± lm pn ^9mp dg nk 



dt dt 



dt dt 



(3.2) 



R 



i k 



g ap R iak p = 9 a0 9kiRL = -9mR\ w + g™g kl R\ 
-Qkl 



lp & 9ip 1 lm pn dcjmp &Qni 

2~dT~dT~2 9 dt 2 ~l 9 9 ~di dt 



+9 m 9ki 
1 d 2 g lk 



r>i , 1 im d 9 mi dg pq 1 lm dg mp dg iq 
n ivi + 4^ 94 ^ 4^ Qt dt 



R 1 dg lp dg ip 1 pn dg kp dg m 1 pq dg ki dg pq 1 pq dgkpdgiq 

2 dt 2 lk 2 9kl dt dt 4 5 9t 0t 4 ff ^ 9t 4 5 dt <9t 



1 a 2 ffifc 



2 5 <9t <9t ■ 



2 9t 2 ** 4^ at at 

A direct computation gives (see, e.g., Fock [4], p. 423; Fisher and Marsden [3], p. 22) 



B - 800 , 1 



1 / ar fc sr fc 

<y + 2 i»« flxi ■ Q x 



(3.3) 



where 



(ft) 



and 



1 —.9 i /37' Le -' 1 — 9 L kl ~ r.9 



kl^O _ 1 kl d 9kl 



klj^l HI "pi " "p 

.9 1 ki — .9 1 fei — 1 



A 



<9t 



l, a g g2 gtf , & / dg a0 \ _l d 2 gij 1 M d 2 gij t q ( dg^p 
2 9 dx«dx? + ij 1 5q/3 ' dx*)~ 2 dt 2 2 9 dx k dx l + ij 1 5q/3 ' 9^ 



Bij ( 3a/3, „ A 



dx 



9 9eatipt ja + - I +5ja1 q/3 3 -Q-r + g lX l af) g g 



1 /% 
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Similar to (3.3), we have 



8M 1 d2 9ij 1 n kl d2 9v | zj f„ 9 5fc' dcta 



ij 2 <9t 2 2 y dx k dx l 'U dt 



and 



,, , dg ki dg kl \ a q ( dg a p 



-<?«it r< + (-i)g kl r° lk r% + g kl : 5p ,rf fc rj z 

1 / <9g»j 1 ki^9ki dgjj k k dg pq p k dg pq 

-2\-dt2 9 -dT+frf^ 9 ' 9 +9ikTrs9 9 ~dxT, 

5fe '2 5 at 2 ff at 5 2 a 2 a +s a* 1 * 1 * 

1 u 9 9ki dg^ I dgjj h 
A 9 dt dt 2 9x fc P<?S 



1 VQ d 9ipd9jq 1 u d 9kl dgjj 
2 9 dt dt 4 9 dt dt 

-a kl a T p V + --hiv k n pq 

9 gp^i^ji + 2 dx ^ Pq 9 

1 I „ T^fc pr qs ®9pq , pfc pr qs ®9pq 



+- 2 {^ rs r^^ + g jk T^g h 

It follows from (3.3) that 

i 'V 2 //,.. , D _ i .v-//,.. i u d 2 9tJ if or or 



2 dt 2 



id* gij i kl d 2 9ij i / sr* ar*\ 

" ~ 2 at 2 2 5 ftrW ' 2 V 1 ' 1 ih- ' <Ar' J 

-ltt („ dgu d 3ki\ , l nPq d9ipdg ]q 1 pq dg i: jdg P q 



i.e., 



da?' dt 2" dt dt 4 a at <9t 



dt 2 +zu v - dt 2 9 dx k dx i + [9ik dxj + 9jk gxi 



(„ d 9ki dg kl \ dg ip dg jq I d gij dg pq 

^{to'Q^'-toJ+a -dT^f-2 9 ~dT~dT 

v 2 9a n M d 2 9ij ( dr^ sr* 
dt 2 9 dx k dx l + \ 9ik dxi 9]k dx- 



+ (gikT^g^ + gjk r k r yr g qsdg^j (JU) 

Similar to the harmonic coordinates in the space-time (see [3]), here we make use of a new kind 
of coordinates on the manifold defined by 

r t g^ri = o. (3.5) 



Such new coordinates are called elliptic coordinates on j$ . 
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Lemma 3.1 Let g^ be a C°° Riemannian metric on the manifold . There is a C°° local 

coordinates transformation <j) : ► , x — ► x around a fixed point p £ ^# such that the 

k 

transformed metric g^ is a C°° Riemannian metric with T (x) = for all x in a neighborhood 
around p £ and any k £ {1, 2, • • • , n}. 

Proof. Consider the elliptic equation for the scalar tp, 



, a h d 2 ip H _,,- dtp 



&X 1 

The coefficients are C°° . Let x 1 (x) be a solution with the condition x l (p) = x % (p) , — — (p) = 5) . 

oxj J 

Then x l is a C°° local coordinates transformation <j> : ./# — ► ^# around p and the transformed 
metric g^ is a C°° Riemannian metric. 

Now the equation Ai[> = is a tensorial (scalar) equation. In the barred coordinates, it becomes 

OXiOXj OXj 

Therefore, satisfies the elliptic condition (3.5). The proof of Lemma 3.1 is complete. % 

By Lemma 3.1, we can choose the elliptic coordinates around a fixed point p £ M and for a 
fixed time t £ R + . After throwing off the bar sign, the geometric hyperbolic flow (3.1) can be 
written as 

d2 9ij _ n kl d2 9ij fr ( n dg kl \ 

~oW~ 9 dxW +Htl { 9kh dxJ)> (3 ' 6) 

where 

ij i 91 * 1 ' ~dx~P J ~ ~ 9 9pg lk jl ~ \ 9 rs9 9 ~dxT 9j rs9 9 ~dx T J ^ ' 

dgki 

are homogenous quadratic with respect to — — and rational with respect to gki with non-zero 

d d 

denominator dct(g) ^ 0. By introducing the new unknowns <fc, hu — gu k — Tr~f: the 

at ax k 

system (3.6) can be reduced to a system of partial differential equations of first order. We now 
consider such a quasi-linear (symmetric hyperbolic) system with ^n(n + l)(n + 2) PDEs of first 
order 



dt ~ hij ' 

ki ®9ijM _ dhjj (3-8) 
dt dx k ' 

dh iQ _ „kl d 9lj,k , fr 

I dt ~ 9 8x1 + ij - 



In the C 2 class, the system (3.8) is equivalent to (3.6). 
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Let u — (gij , gij.k, hij) T be the — n{n + + 2)-dimensional unknown vector function. The 
coefficient matrices A a ,A J ,B are given by 



A°(u) = A (gi j ,g ij , k ,hi j ) = 



( I 











o\ 





9 n I 


9 12 I ■ 


■ 9 ln I 








9 21 I 


9 22 I ■ 


■ 9 2n I 





o 


a nl I 
y 1 


a n2 I ■ 
y x 


■ a nn I 
y 1 


o 


V 











1 J 


o 


o 








o 














9 3l I 














9 J2 I 














9 3n I 





9 lj I 


ir'i ■■ 


• 9 nj I 






A J (u) = A J (g ku gki,p,h k i) 



I 

where is the ( —n(n + 1) I x ( —n(n + 1) I zero matrix, / is the ( —n(n + 1) I x ( —n(n + 1) 



identity matrix, 



B(u) = B(g i: j,g. ljjP , h i3 ) 



( , 





in which is the —n 2 (n + 1 )-dimensional zero vector. 
We observe that the symmetric hyperbolic system 



(3.9) 



is nothing but the system (3.8). So far, we have reduced the equation of the hyperbolic geometric 
flow (3.1) to the symmetric hyperbolic system (3.9), which are equivalent to each other in the C 2 
class. Then, by the theory of the symmetric hyperbolic system, the smooth solution to (3.1) exists 
uniquely for a short time (see [3]). Thus, the proof of Theorem 1.1 is completed. 

Remark 3.1 The elliptic coordinates can also be used to prove the short-time existence for the 
Ricci flow. 

More generally, motivated by general Einstein equations and the rich theory of hyperbolic 
equations, we may also consider the following field equations with the energy-momentum tensor 
Tn under certain conditions: 



a gij 



2i? ? ; 



9 9 \ _ r 



(3.10) 
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where are certain smooth functions on ^# which may also depend on t, #y are some given 
smooth functions of the Riemannian metric g and its first order derivative with respect to t, and 
K is a parameter. Similar results can be obtained. 

4 Nonlinear stability for hyperbolic geometric flow 

In this section we investigate the nonlinear stability of the hyperbolic geometric flow defined on 
the Euclidean space with the dimension larger than 4. 

We now state the definition of nonlinear stability of the hyperbolic geometric flow (1.1). Let 
^ be an n-dimensional complete Riemannian manifold. Given symmetric tensors <A and g}j on 
we consider the following initial value problem 
( d 2 

— gij(x,t) = -2R i:j (x,t), 

9ij{x,0) = g i3 {x) +eg% j (x), -gf(x,0) = egj^x), 
where e > is a small parameter. 

Definition 4.1 The Ricci flat Riemannian metric ~g~ij{x) possesses the (locally) nonlinear stability 
with respect to (gfj,g}j), if there exists a positive constant Eq = £o(gfj,gjj) such that, for any 
s € (0,£o], the initial value problem (4-1) has a unique (local) smooth solution gij(x,t); 

~g~ij(x) is said to be (locally) nonlinearly stable, if it possesses the (locally) nonlinear stability 
with respect to arbitrary symmetric tensors gfj(x) and g}j(x) with compact support. 

In what follows, we consider the nonlinear stability of the flat metric of the Euclidean space 
R™ with the dimension n > 5. We have 

Theorem 4.1 The flat metric gij — 5{j of the Euclidean space R™ with n > 5 is nonlinearly stable. 

Remark 4.1 Theorem 4-1 gives the nonlinear stability of the hyperbolic geometric flow on the 
Euclidean space with dimension larger than 4- The situation for the 3-, 4~dimensional Euclidean 
spaces is very different, and will be studied in the sequel by using null conditions. This is similar to 
the proofs of the Poincare conjecture in topology: the proofs for the three and four dimensional case 
and n > 5 dimensional case are very different (see, for example, JfJjj and Jj]/). This motivates us to 
understand the possibility of using hyperbolic geometric flow to understand Poincare conjecture. 

Proof of Theorem 4.1. Define a 2-tensor h in the following way 

gij(x,t) = % + hij(x,t). 
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Let 5 l i be the inverse of Sij . Then, for small h 



where h l i = S lk Si l hki and O tJ (h 2 ) vanishes to second order at h = 0. 

As in Section 3, we choose the above elliptic coordinates {x 1 } around the origin in W 1 . Then 
the initial value problem (4.1) can be written as 



(*,*)= (*« + H-) *^ + H« (S kl + h kl (x, t) 0hl - 1 ' 1 



d 



hij(x,0) = egfj(x), — %(x,0) = egjj(x), 



where 



dxP 



(4.2) 



Hij S ki + hki(x,t), 



dh k i{x,t) 
dxP 



Hij gki(x,t) 



dg k i(x,t) 
dxP 

,9g Pg 



dxJ 



dx i _ 
dh 



-2 (S kl + IT*) (6 pq + h pq ) - (5 pa + H pa ) (S qb + H qb ) , } . 
- (S ik + H ik ) (6 pr + H pr ) (6 qs + H qs ) 1 (5 ka + H ka ) (^f 



1 / s-i It 

- {8 jk + Hj k ) (S pr + H pr ) (S qs + H qs )- (S ka + H ka )' 



dh a k dhik 



dx % 



-\s ik 8 pr 5 qs 5 ka {^- s 
dh„. 



dh a 



dx a 
dh rs \ dh 



\ dx s 
dh b: j dhu 
dx l dxi 



S jk 5 pr S qs S- 



r ^qs %ka 

ab (dh, 



1 

~2 

-h kl s c , 

2 V dx 

--§p r s qs ( dhir 

2 V dx s 

2 V dx s 



dx r 
dh a 



dx a J dxi 
dh rs \ dh pq 



+ 



dx s dx r dx a J dx 1 
dh ak dh ik 



O \h 



dx 1 
dh ls 
dx r 

dh is 



dx" 
dh rt 
dx % 
dh r 



dh bj dhu 
+ 



dh 



dx 1 dxi Q x b J 



dh 



pq 



dx r dxi 
Furthermore, (4.2) can be written as 



dxi 
dh pq 
dx % 



0[\h 



i ki 



dh k 



dxP 



dh a k 


dhik^ 


/ dh bj 


dx 1 


dx a j 


\ dx 1 


dh as 


dh rs \ 


dh pq 


dx r 


dx a J 


dxi 


dh as 


dh rs ^ 


\ dhpq 


dx r 


dx a , 


1 dx 1 


dh 3 i \ 






dx b ) 







+ 



dh 



hi 



dxP 



dhu _ dhji 
dxi dx h 



' ' hij(x,t) =J2~ — 



dt 2 



fe=i 



dx k dx k ^ 



hkl, 



dh k i d'' 



hu \ 



dxP ' dxPdxi / ' 



h tj (x, 0) = sgfj (x) , — h^ (x, 0) = eg\j (x) , 



(4.3) 
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where 



dh ki d 2 h ki 



HlJ \ kkh dxP ' dxPdxi J 



\ = 



H 



d 2 h 
kl 67 'Hj 

dx k x l 



Ha S k i + h 



Oh 



kl 



'77 



— 5 kl S a 
2 



2 



2 





dh ak 


dx k 


dx l 


dhi r 


dhi S 


dx s 


dx r 


dhj r 


dh js 


dx s 


dx r 



' dxP 
dh ik 
dx a 
dh rs 
dx l 
dh rs 



(dh 



bj 



dh 



\ dx l 
dh pq 
dxi 
dh pq 
dxi I dx l 



dh h 



dxi dx b 



k:l 



d 2 h i:j 
dx k dx l 



O 



\\h u \ + 


dh u 


+ 


d 2 h u 




dxP 




dx p dx q 


) 



Let 

- / dh u d 2 h H \ . , . _ 
X= \ hkh ^^d^) iP^k,l = l,2.. 

The nonlinear term 

^ {X) = ^{ hkl ^^jFxr 

in (4.3) is smooth in a neighborhood about A = and satisfies 



,n). 



H i3 {\) = 0[\\\ z ) = 1,2,- ■■ ,n). 

By the well-known global existence results for the nonlinear wave equation (e.g., see [2], [7]> [S], 
|10j). there exists a unique global smooth solution (hij (x, t)) for the Cauchy problem (4.3) or (4.2). 
Thus, the proof of Theorem 4.1 is complete. 



5 Wave character of the curvatures 

The hyperbolic geometric flow is a system of hyperbolic evolution equations on the metrics. The 
evolution of the metrics implies a system of nonlinear wave equations for the Riemannian curvature 
tensor Rijki, the Ricci curvature tensor Rij and the scalar curvature R which we will derive. 

Let ^# be an n-dimensional complete manifold. We consider the hyperbolic geometric flow on 
. // . that is, 

d 2 

g^9ij(x,t) = -2Rij(x,t). (5.1) 

We now want to find the evolution equations for the Riemannian curvature tensor Rijki, the Ricci 
curvature tensor Raa and the scalar curvature R. 
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Direct computations yield 



dt 2 jl 



p h 

"■ijl 

d 2 

° p h 

Q^2 ijl 

d^ 

Q£2 Rijkl 



I hm 

2 y 



' dg mj 

y 8x l 



dg m i _ dgji \ 
dxi dx m ) '' 
d 2 g m i _ d 2 g jt 



\ n h m ( 9 2 g m j d 2 g m i _ 

^dx l dt dxidt dx m dt 



2 



1 d 2 g hm 



dg mj ogm, 

\ dx l dxi 
( 9 ( d 2 g mj 
\dx l V dt 2 
fir h 

il , pli rP _ ph - 

<).,- dxi ^ v H 3V 
,9a; 1 \dt 2 jl J dxi \dt 2 



2 dt 2 

+\g hm 



d f d 2 g mi s 
dxi ' 



i hm ( dgmj dg m i _ dgji \ 
dt \ dx l dxi dx m ) ' 

dx m dt y 



dgmi dg jt \ 1 dg hm ( d 2 g mj d 2 g mi 
dx m ) '2 dt \dx l dt dxidt 



dt 2 



d 

dx m 



dt 2 



+ dt 2 



~\h j-\h 



ph pp 

jp il 



dt 2 

9hk 



1 dghk 
+I ~dT 







_dx l 


\dt 2 ^ 1 ) 



d 

dx l 



(■pi- , 
V x ip L jl 

d 2 r>h , D h d 2 ghk dghk d uh 
dxi \dt 2 il ) 

L (In 



' dt 2 



'ph pp 
1 ip L jl 



■' 1 dxi \dt 



' + I( r W 



We choose the normal coordinates around a fixed point p on such that 

r&(p) = o, 

or, equivalently, 

dx k 



^i) 

-?i P n) 



, T>h d 2 ghk 
+ K iji-Qtr- 



1G 



Then we have 



Q-f.2 "^jkl 



9hk 



dx 1 



1 d 2 g hm 

2 dt 2 



dg r , 



dgmi dg 3 i 



+ghk 



'9hk 



a 

d 
dxi 

d 

9hk dxj 

+2ghk 
dghk 



1 



7-9 



hm 



dx l dxi 
d fd 2 



2" \dx l 
1 dg hm (dg 

1 



Qmj 

dt 2 

dg 



dx r > 
d 

dxi 



ldg hm f d 2 g mj 
' 2 dt { dx l dt 



d 2 g m i d 2 gji 



dx^dt dx m dt 



d 2 g m i 
dt 2 



da 



_ dgu 
dx 1 dx m 



ml 



2 9 



hm 



_d_ 

dx 1 



d 2 g m 
dt 2 



+ ■ 



d ( d 2 g m i 



dx 1 \ dt 2 



d 

' dx m 

1 dg hm 

2 dt 

d 



'9ji 



dx r > 



dt 2 

d 2 a ■ 
dx l dt 

d 2 g l i 
dt 2 



d 2 g m i _ d 2 g a 
dx l dt dx m dt 



dt ip dt i l dt ]P dt il 



+2 



dt 



1 d 

2dx l 



dg hm dg, 



r, d 9hk 1 hm d 
dt 2 y dx 1 



, dg h k 1 

' dt 2' 



Jim 



1 Q 2 g hm 

ghk- 



2 

+ghk 
+ghk 



dt V dx 1 

d ( dg m j 
dx 1 \ dt 

d_ ( d_ ( dg mi 
dxi ydx 1 \ dt 

dg, 



dg m i 
dxi 

d 

dxi 



d 9jl\\ 
dx m ) ) 



1 d 

2 dx^ 



dg hm (dg, 
~dt 



dx 1 



mi ^ dg m i 



dx 1 



dgu 

Q x m 



dg m i\ 


d 


( dg 3 i 


dt ) 


dx m 


\ dt 


dg m i\ 


d 


f dgu 


dt 


dx m 


{ dt 



+ R: 



d 2 g h k 



_d_ 

dx 

(dg 



dx 1 



dt 2 

02 ghm 

dx l dt \ dx'dt ' dxidt 

dg hm 
dt 

d 2 



%! 

dxi Q x m 

dgji 



d 
dxi 



Qghm / Q g 



dt \ dx 1 
(dg 



+ 



ijl g t 2 

dg m i dgu 



dx 1 



d_ (dg 

dxi ydxigt ' dxidt dx m dt 



mj $Qml 



dx m dt 
dgji 



q2 ghm 

9hk dxidt \ dx'dt dx* dt ~ dx m dt 



'm.i $9 ml 



dx" 
dgu 



d ( d 2 g mi d 2 g m i 



dx l dx l 

d 2 



d 2 gkj 
dt 2 

d 2 gki 
dxidx 1 V dt 2 



d 2 



dx % dxi 

d 2 

dxidx 1 
£> . d 



d 2 g k i 
dt 2 

d 2 g k i 
dt 2 



dxi y dx l dt 

d 2 f d 2 9]l 
dt 2 

d 2 gu 



d 2 gu 

dx'dt dx m dt 



dx l dx k 

d 2 



dxidx k \ dt 2 



+2ghk I dt T ^ dt vP > 1 dt r ^dt Til 



h m 



+ 



dghk dg 
dt dt 



d 

dx 1 



dg 



dg m i dgji \ 



dx 1 



I 9 9hk hm 

dt y 



d_ (dg 



dx* \dx l dt + dxidt dx m dt) dxi \dx l dt^ dx l dt dx m dt 



dxi 
dg m i 



dx m J 
dgji 



d 

dxi 



dg 



mi , dg m i 



dx 1 dx 1 



dgu 
dx m 



d ( dg mi : dg mi 



+ 



dgu 



+R: 



h d 2 g h k 
ijl dt 2 ■ 



(5.2) 
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Noting g hm g ml 



Sj 1 , we get 

Q g hm ^ 

dt = 

d 2 g hm 

dx k dt 

Ql g hm 
dt 2 

Thus, it follows from (5.2) that 

d \ ( l aP md 2 gu P 



hp mq d 9vi 
99 dt ' 

hy mq d 9pq 

9 9 dx k dt 

hy mQ d2 9pq 

9 9 dt 2 



+ 2g hp g rq g sm - 



dg pq dg r 



dt dt 



2 

_d_ 

dx 



dt 2 

dg m j 
dx 1 



+ g rq g pm 



rp 



+ 



dg, 



ml 



dxi 



,d 2 



: vm - 2 9kp f d 2 g m j 
dx l dt \ dx l dt 



d 2 



+2ghk 



dx l dx l 

d 2 

dxidx 1 
d 



dt ip 



d 2 gkj 
dt 2 

d 2 gti 
dt 2 

d_ r p 
' dt * 



dgkg dg. 
dt dt 

dgji \ 
dx m ) 

d 2 g m i _ d 2 g jt 
dxidt 

d 2 



d 

dxi 



dg 



mi ^ dgrfil 

dx 1 dx 1 



dgu 
dx m 



+ 



dx l dxi 

d 2 

dxidx* 
d 



dx m dt 

d 2 g k i 
dt 2 

d 2 g k i 



, pin 



d 2 gk P ( d 2 g mi d 2 g m i d 2 g a 



dxidt V dx l dt dx l dt dx m dt 



d 2 



dt 2 



dt ]P 



d_ r p 

dt 11 



R 



—g hp g mq 



d 2 gk P 

dt 2 
d 2 gk P 



dg n k dg pq 
dt dt 

ArP - 

dx* jl 



d 

dx 1 

d 
dxi 



1 .;; 



dg m j dg, 
dx 1 

+ 2g rq 



ml 



dxi dx 
dgkg dg rp 



dx l dx k 

d 2 

dxidx k 

l d 2 g hk 
ijl g t 2 

) 



d 2 9pi 
dt 2 

d 2 gzi 
dt 2 



_ dg jt \ d_ 

1 dxi 



dg m i , dg 



ml 



dx 1 dx 1 



dgu 
dx m 



dx l dt 



+ 



d 2 



+ighk 



dx l dx l 

d 2 

dxidx 1 
d 



dt 



, ( d 2 9mj 
V dx l dt 

(-2R kj ) + 

(-2Rik) + 

d_ r p 

dt * 



d 2 g m i 
dx^dt 

d 2 



dt 
fi 9ji 



dt 



_d_ r p 

dx* 11 



JLvp 

dxi « 



dx m dt 



9kp ^ gP ni 



( d 2 g mi d 2 g ml d 2 g a 



dx l dxi 

d 2 



dx % dxi 

*I* ■ 

dt 1P 



- (-2Rki) - 
(-2Rki) - 



dxidt 
d 2 



\dx l dt dx*dt dx m dt 



dx l dx k 

d 2 

dxidx k 



(-2%) 
(-2Ril) 



d_ T P 

dt a 



dghk dg. 



dt 



d 2 



dx l dx l 

d 2 



(-2R kj ) + 



d 2 



{-2R ki ) + 



dx % dxi 

d 2 



dxidx 1 

d 2 gk P ( d 2 g m j 



+ 



dx % dxi 
d 2 g m i 



(~2R kl ) - 
(-2R k i) 



dt 

d 2 

dx l dx k 

d 2 



Elghp . 2 



d 



n, 



dxi jl 



dxj %l 



R 



d 2 ghk 
ijl dt 2 



d 2 9ji 



dx % dt \ dx l dt dxidt dx m dt 



dxidx k 
d 2 



(-2%) 
i-2Ra) 



+ g 



pm _ 



( d 2 g mi d 2 g mi d 2 g a 



+2ghk 



—v h 

dt ip 



d T P 

dt i l 



dt 3P 



9_ t p 

dt a 



dxidt \ dx'dt dx*dt dx m dt 



(5.3) 



On the one hand, we have 

d 2 



dx l dx l 



R 



jk 



ViViR jk + v,rf fc • R jp + v^. 



Rkp- 
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Then 



d 2 



dx l dx l 

d 2 



(-2R k j) 



d 2 



{-2Rki) + 



dx % dxi 

d 2 



(-2R k i) - 
(~2Rki) 



d 2 



dx l dx k 

d 2 



(-2%) 
(-2R U ) 



dxidx 1 dx l dx'i dxidx k 

— VjVii?^- - ViT p lk R 0P - ViTfjR kp — Vi\7jR k i — ViT^Rip - V^T^i?^ 
+V l V k R J i + ViT p kj R lp + V t T p kl R 3P + V l W i R kt + V 3 T p kl R ip + V 3 T p u R kp 
+V 3 V l R kl + VjT p ik R lp + V 3 YlR kp - V 3 W k R a - V^.i^ - V ,T p kl R pi 
—ViViR k j — ViVjRki + ViVkRji + VjV;i?fci + VjViRki — VjV k Ru 

+R iP (Vj-rffc - v,r^) + R jp (-v,rf fe + v.r* ) 

+R kp (-v^ - + v,-rf. + v,-rs) + R lp (-v t r P k + v^. + v,rf fc - v,r^) 

—Riji P g pq R q k — Rijk P g pq Rqi + Rk P 2-R^/) 

-ViViRkj + ViVfciZj-j + Vj-VjiZfci - Vj-Vfci2« - 5 P9 {Riji P R kq + Rijk P Ri q ) ■ 



On the other hand, we have 
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d 2 fffc P / d 2 g mj d 2 g mt _ d 2 g 3 i 



dx*dt V dx l dt + dxidt dx m dt 

^V h —V p - —V h —V p 

dt ip dt * dt jp dt a 



+ g p 



d 2 gk P f d 2 g mi d 2 g m i _ d 2 g lt 



d 2 g kp ( d 2 g mj d 2 g m i _ 
dx'dt { dx l dt dxidt dx m dt 



9ji 



+ g p 



dxidt \ dx l dt dx l dt dx m dt 



, d 2 g kp / d 2 g ml d 2 g ml _ dP_ 



dxidt \dx l dt dx l dt dx m dt 



+ -g 
2 y 



1 pm ( d 2 g k i d 2 g kp _ d 2 g ip \ ( d 2 g mj d 2 g mi _ d 2 g jt 



\ dxPdt dx l dt dx k dt J \ dx l dt dxidt dx m dt 



2 y 



( d 2 gk 3 



+ 



d' 



gk P 



d 2 



gjp \ ( d 2 g ml d 2 g m i d 2 g a 



pm 



\dxPdt dxidt dx k dt J \dx l dt dx^dt dx m dt 
1 / d 2 g mj d 2 g„u d 2 g jt \ ( d 2 g ki d 2 g ip d 2 g kp 



2\dx l dt dxidt dx m dt J \dxPdt dx k dt dx l dt 



1 / d 2 g ml d 2 g mt _ d 2 g a \ ( d 2 g kj _ d 2 g JP 



2\dx l dt dx'dt dx m dt ) \dx?dt dx k dt dxidt 



gk P 



)} 



= 2 5p 



1 „ P r 9 ( d 9rj , dgn dgji \ 1 nQS d ( dg kl dg ls dg ks 



T-9 



dt V dx l dxi dx r 



qs_ 



dt \ dx s dx k dx l 



- a w!L ( ^ili + QtL - . - g i s — ( ^M. _ 



2" dt \ dx l 
2g pq 



dg 



ks 



dx % dx r ) 2 dt \ dx s dx k dxi 



£Lr p —T q - —V p —V q 

dt ll dt 3k dt i l dt lk 
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Therefore, it follows from (5.3), (5.4) and (5.5) that 
d 2 

-gpRijkl = — VjVi-Rfcj + ViVkRjl + VjVlRki — VjVkRil — 9 P<1 {RijqlRkp + RijkqRkp) 



dt a dt ]k dt jl dt u 
Similar to Hamilton [6], we have 

Theorem 5.1 Under the hyperbolic geometric flow (5.1), the Riemannian curvature tensor Rijki 
satisfies the evolution equation 
d 2 

-QpRijki = ARijki + 2 (Bijki — Bijik — Bujk + Bikji) 

—g P<1 (RpjklRqi + RipklRqj + RijplRqk + RijkpRql) (5.6) 

^Lr p — v q -— v p — v q 
dt a dt i k dt i l dt ik 

where Bjjjy — g pr g qs R p i q j Rrksi and A is the Laplacian with respect to the evolving metric. 



Remark 5.1 In Theorem 5.1 and Theorem 5.2 below, the term 2q n „ | — T p , ■ — T q , — — T p , ■ — T q , 

' zpq I a. il q. ift q. jl o, ik 



Ar" —T q ~—Y p — 
dt a dt jk dt 11 dt^ i 

can be written in the covariant form. For the sake of simplicity, we omit it. 

For the Ricci curvature tensor, we have 

d 2 d 2 

dp Rlk = dt? i^kig 3 ) 

d 2 d d d 2 g^ 

— r,3 l ® R n„ivjQ ®9P1 ® P „:ipjq d 2 9pg p i r, „jp „rq „sl ®9pq ^drs D 

- 9 W R llk i~2g g —-^-g g —R ljkl +2g g g ——R l]kl . 
Thus, we obtain 

Theorem 5.2 Under the hyperbolic geometric flow (5.1), the Ricci curvature tensor satisfies 

d 2 



dt 2 



Rik — Ai?jfe + 2g pr g qs R piq kRrs — 2g pq R pi R q k 

+V*.(|lS|ll-|ls|ll) (5.7) 



For the scalar curvature, we have 

d 2 d 2 

Qf2 R = -^(9' ^k) 



ik d 2 R , 9 d d k d 2 g ik 

k —R. r, n ip n kq ®9pq dRik p / „ip „kq 3 2 9pq , r, jp n rq n sk ®9pq ^9rs 



n lK T3 Q„tp „Kq - zpq '» i r> ip n kq "JM , ip rq 

9 g^Rik-2g g — — + g -^ + 2g g g ^ ^ 
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On the other hand, 

2a ik a jl a (— T p — T q - —T p —T q 
9 9 9vq \dt il dt » dt ^dt ik 

= f*VV'V r (^)v.(2gL) - i 5 "V r (9 ife ^)V s (^^) 

+2A^ r ( 9 ' t f)V I (^)-//fV r (|)V I (|i)-2 9 Vs"V I (^)V J (|). 
Then, we get 

Theorem 5.3 Under the hyperbolic geometric flow (5.1), the scalar curvature satisfies 
d 2 

= Ai? + 2|Ric| 2 

at 

+2 5 '-V'V r (5 <fc ^)V I (^ ) - ffVV a V r (^)V,(^) (5-8) 
-2g gig V^—)V 3 ( — )-2g fPg*—-^ 



at dt ' ™» * * dt dt ■ 

Theorems 5.1-5.3 show that the curvatures of the hyperbolic geometric flow possess the wave 
character. We will apply techniques from hyperbolic equations to the above wave equations of 
curvatures to derive various geometric results. 



6 Discussions 

The hyperbolic geometric flow describes the wave character of the metrics and curvatures of mani- 
folds. Many hyperbolic systems in nature provide natural singular sets, the typical example is the 
Einstein equations in general relativity which form a hyperbolic system with a well-posed Cauchy 
problem. If one starts with smooth initial data, one may end up with a singular space-time. One 
of the most challenging problems is to describe the kind of natural singularity. The famous cosmic 
censorship conjecture due to Penrose is an attempt to describe such singularities (see [12 ). In 
Kong and Liu [TT] , we construct some exact solutions of the hyperbolic geometric flow, these solu- 
tions possess the singularities which are nothing but those described by Penrose's conjecture. By 
these examples, we believe that the hyperbolic geometric flow is a very natural and powerful tool 
to understand the singularities in the nature, in particular, the singularity described by Penrose 
cosmic censorship conjecture. 

The Einstein equations play an essential role in general relativity. Consider a space-time with 
Lorentzian metric 

ds 2 = g^dx^dx 1 ' (/x, v = 0, 1, 2 • • • ,n). (6-1) 
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The vacuum Einstein equations read 

Gfiv = 0, (6.2) 
where G^ v is the Einstein tensor. We now consider the following metric with orthogonal time-axis 

ds 2 = -dt 2 + gij (x, t)dx i dx j . (6.3) 

Substituting (6.3) into (6.2), we can obtain the equations satisfied by the metric 

dt 2 - zn v 2 9 dt dt +9 dt dt • [b ' 

Neglecting the lower order terms gives the hyperbolic geometric flow (1.1). Therefore, in this sense, 
the hyperbolic geometric flow can be viewed as the leading terms in the vacuum Einstein equations 
with respect to the metric (6.3). Since the hyperbolic geometric flow only contains the main terms 
in the Einstein equations, it not only becomes simpler and more symmetric, but also possesses 
rich and beautiful geometric properties. In particular, in mathematics, its Cauchy problem is well- 
posed and easier to handle some fundamental problems such as the global existence and formation of 
singularities; on the other hand, it can be applied to re-understand the singularity of the universe 
and other important problems in physics and cosmology (see |16j). We also believe that there 
should be some relations between the solutions of the Einstein equations and the corresponding 
hyperbolic geometric flows. On the other hand, from the above discussions we have seen that the 
hyperbolic geometric flow also possesses many beautiful features similar to those of the Ricci flow, 
and some of the techniques in the study of the Ricci flow can be directly used to understand the 
hyperbolic geometric flow. The deep study on the hyperbolic geometric flow may open a new way 
to understand the complicated Einstein equations. 

It is well known, in general relativity there is a constraint system of equations involving an 
asymptotically flat metric tensor and another symmetric tensor. There are four constraint equa- 
tions and it is therefore over-determined. Unlike this, since the time axis is orthogonal to other 
space axes, the hyperbolic geometric flow does not need to satisfy any additional constraint. More 
precisely, for the Cauchy problem of the hyperbolic geometric flow, in order to determine the solu- 
tion we need two initial conditions: one is the metric flow itself gij(x,0), another is its derivative 
^jr-(x, 0), since the time axis is orthogonal to other axes, these initial data do not need to satisfy 
any additional constraint, and therefore it is a determined system. This is another main new 
feature of the hyperbolic geometric flow. 

Many mathematicians, for example Shatah et al [13j-|15j. have investigated the Cauchy problem 
for some geometric wave equations. The model at hand is the harmonic map problem, which is 
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the study of maps from the Minkowski space-time into complete Riemannian manifolds. This kind 
of geometric wave equations is a system of partial differential equations of second order, which is 
the Euler-Lagrange equations of the action integral of the harmonic map. It satisfies certain linear 
matching condition, and then under suitable assumptions, has a unique small smooth solution 
for all time, and possesses some interesting (decay, energy and regularity) estimates. On the 
other hand, the hyperbolic geometric flow is determined by the Ricci curvatures of a family of 
Riemannian metrics on the manifold under consideration. That is to say, the hyperbolic geometric 
flow possesses itself intrinsic geometric structure and can be used to describe the wave character 
of metrics and curvatures. This is essentially different from the above harmonic map problem. 

As well-known, one can understand the heat kernel from the kernel of wave equation. This 
indicates that we should be able to derive various information of the Ricci flow from that of the 
hyperbolic geometric flow. Therefore it is also interesting to understand the relations between the 
hyperbolic geometric flow and the Ricci flow, the singularities of its solutions and its relation with 
the geometrization theorem. This will be another interesting topic in the sequel. 
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